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We present a formalism to study the metric perturbations of the Schwarzschild spacetime. The for- 
malism is gauge invariant, and it is also covariant under two-dimensional coordinate transformations 
that leave the angular coordinates unchanged. The formalism is applied to the typical problem of 
calculating the gravitational waves produced by material sources moving in the Schwarzschild space- 
time. We examine the radiation escaping to future null infinity as well as the radiation crossing 
the event horizon. The waveforms, the energy radiated, and the angular-momentum radiated can 
all be expressed in terms of two gauge-invariant scalar functions that satisfy one-dimensional wave 
equations. The first is the Zerilli-Moncrief function, which satisfies the Zerilli equation, and which 
represents the even-parity sector of the perturbation. The second is the Cunningham-Price- Moncrief 
function, which satisfies the Regge- Wheeler equation, and which represents the odd-parity sector of 
the perturbation. The covariant forms of these wave equations are presented here, complete with 
covariant source terms that are derived from the stress-energy tensor of the matter responsible for 
the perturbation. Our presentation of the formalism is concluded with a separate examination of 
the monopole and dipole components of the metric perturbation. 



I. INTRODUCTION 

Metric perturbations of the Schwarzschild spacetime 
have been studied for a long time, starting with the pio- 
neering work of Regge and Wheeler pj , Vishveshwara , 
and Zerilli 0. The theory was summarized in an influ- 
ential monograph by Chandrasekhar |4| , and it has been 
applied to many different physical situations (see, for ex- 
ample, Chapter 4 of the book by Frolov and Novikov 
In particular, a useful application has been the compu- 
tation of gravitational waves produced by a point par- 
ticle moving in the field of a Schwarzschild black hole 
(see Ref . , the review article of Ref . , and references 
therein). Another has been the simulation of a collision 
of two black holes in a "close-limit" approximation (see 
the review article of Ref. @ and references therein) . 

Traditionally the perturbation formalism is developed 
in the standard Schwarzschild coordinates (t, r, 9, </>), and 
in a standard choice of gauge known as the "Regge- 
Wheeler gauge." The tradition also makes use of Fourier- 
transform techniques and presents the perturbation equa- 
tions in the frequency domain instead of the time do- 
main. Moncrief [!j was the first to present the formalism 
in a gauge-invariant package, recognizing the practical 
advantages that gauge invariance provides: While the 
Regge- Wheeler gauge is useful for many purposes, it is 
not useful for others, and the power to switch from one 
gauge to another within a gauge-invariant framework is 
often required. Another refinement of the formalism was 
produced by Gerlach and Sengupta 10], who presented 
the gauge-invariant perturbation equations in an arbi- 
trary coordinate system, thus liberating the formalism 
from the usual Schwarzschild coordinates and their poor 
behavior at the event horizon. 

This program to translate the traditional perturbation 
formalism into a covariant, gauge-invariant language was 
recently revived by Sarbach and Tiglio |Tlj . and with 



this paper we hope to deliver its final chapter. We aim 
to present the formalism of metric perturbations of the 
Schwarzschild spacetime in its most mature and prac- 
tical incarnation yet. Our formalism is covariant and 
gauge invariant, and it goes beyond the work of Ger- 
lach and Sengupta, and of Sarbach and Tiglio, by the 
development of covariant master equations for the even- 
parity and odd-parity sectors of the theory, complete with 
covariant source terms that are derived from the stress- 
energy tensor of the matter which is responsible for the 
perturbation. 

We have in mind a typical application of the per- 
turbation formalism, the calculation of gravitational 
waves produced by material sources moving in the 
Schwarzschild spacetime. We are interested in the ra- 
diation that escapes to future null infinity and manifests 
itself as waveforms h + and h x that are directly observ- 
able to gravitational-wave detectors, and we are inter- 
ested in the radiation that crosses the black-hole horizon. 
This application illustrates well the need for a covariant 
and gauge-invariant formalism, as the two types of radia- 
tion admit different descriptions. The radiation at future 
null infinity is best described by casting the perturba- 
tion in an outgoing radiation gauge and expressing it in 
a retarded coordinate system (u, r, 9, </>) related to the 
usual Schwarzschild coordinates by the transformation 
u = t-r- 2Mln(r/2M - 1). The radiation at the hori- 
zon, on the other hand, is best described by adopting an 
incoming radiation gauge and expressing the perturba- 
tion in an advanced coordinate system (v, r, 9, 4>) related 
to the usual Schwarzschild coordinates by the transfor- 
mation v = t + r + 2Mln(r/2M — 1). Our formalism 
permits the use of any coordinate system x a — (a; , a; 1 ) 
that can be obtained from the usual Schwarzschild coor- 
dinates (i,r); we do not, however, consider transforma- 
tions of the angular coordinates 9 A = (9, <j>). 

In the formalism developed in this paper, the radiation 
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at future null infinity and at the horizon are described 
in terms of scalar and gauge-invariant master functions, 
^evenC 2 -") and ^f l ^ d (x a ), which can be computed from 
the metric perturbations. These functions are labeled by 
spherical-harmonic indices I and m, and by their behavior 
under a parity transformation. The function ^> l ™ en (x a ) 
is constructed from the even-parity perturbations, and 
it is equal to the gauge-invariant function that was first 
introduced by Moncrief Hk it is a close cousin to Zerilli's 
original master function |3|, and it satisfies a covariant 
version of Zerilli's differential equation. The complete 
source term for this equation is presented for the first 
time in this paper. The function ^> ™ d (x a ), on the other 
hand, is constructed from the odd-parity perturbations, 
and it is equal to the gauge-invariant function that was 
first introduced by Cunningham, Price, and Moncrief [l2T ] 
and recently revived by Jhingan and Tanaka |l3j| ; it is es- 
sentially the time integral of the original Regge- Wheeler 
master function 1| , and it satisfies a covariant version of 
the Regge- Wheeler equation. The complete source term 
for this equation is presented for the first time in this 
paper. 

The paper is organized as follows. In Sec. [n] we give a 
covariant description of the Schwarzschild spacetime and 
specify our notations and conventions. In Sec. I11II we in- 
troduce the scalar, vector, and tensor spherical harmon- 
ics that are used in the decomposition of the metric per- 
turbation. In Sec. IIVI we examine the even-parity sector 
of the perturbation, introduce the Zerilli- Moncrief master 
function, and derive the one-dimensional wave equation 
that it satisfies. In Sec. [V] we examine the odd-parity 
sector of the perturbation, introduce the Cunningham- 
Price-Moncrief master function, and derive the one- 
dimensional wave equation that it satisfies. In Sec. IVT1 
we describe the behavior of the perturbations near fu- 
ture null infinity, construct the radiative part of the per- 
turbation field, extract the waveforms h + and h x , and 
compute the rates at which the radiation carries away en- 
ergy and angular momentum. In Sec. lVIll we describe the 
behavior of the perturbations near the black hole's event 
horizon, and calculate the rates at which they transfer 
energy and angular momentum to the black hole. In 
Sec. IVIlTl we conclude our presentation and examine the 
nonradiating pieces of the metric perturbation, those as- 
sociated with the low multipoles / = and 1 = 1. 

Various technical details are relegated to the Appen- 
dices. In Appendix [S] we expand our discussion of vec- 
torial and tensorial spherical harmonics. In Appendix 
IbI we present the perturbed Ricci tensor for a general 
spherically-symmetric background spacetime. In Appen- 
dices O |D| and [E] we list the perturbation equations 
in the commonly used coordinate systems x a = (t,r), 
x a = (it, r), and x a — (v,r), respectively. 

Because the topic of metric perturbations of the 
Schwarzschild spacetime is so venerable, we will allow 
ourselves in this paper to simply state our results and 
omit most lengthy derivations that lead to those results. 
We hope, however, that the path to the results will always 



be clearly delineated. We refer the reader to the litera- 
ture reviewed in this Introduction for additional details; 
another repository of relevant derivations is Martel's PhD 
dissertation 01 . Throughout this paper we adopt the 
sign conventions of Misner, Thorne, and Wheeler [Tlj . 
and we set c = G — 1 . 



II. SCHWARZSCHILD SPACETIME 

The Schwarzschild metric is expressed as 

ds 2 = g ab dx a dx b + r 2 n AB dO A dO B , (2.1) 

in a form that is covariant under two-dimensional co- 
ordinate transformations x a — > x' a . The coordinates x a 
span the submanifold M 2 of the Schwarzschild spacetime 
- the "(t, r) plane" — and lower-case Latin indices a, 
b, c, etc. run over the values and 1. The coordinates 
9 A = (9, <fi) span the two-spheres x a = constant, and 
upper-case Latin indices A, B, C, etc. run over the values 
2 and 3. The full spacetime manifold is M = M 2 x S 2 . 
The two-dimensional tensor g a b and the scalar r are func- 
tions of the coordinates x a , and £Iab = diag(l, sin 2 9) is 
the metric on the unit two-sphere. 

We shall use three different coordinate systems x a in 
the applications of the perturbation formalism to be pre- 
sented below. The first is (t, r), the usual Schwarzschild 
coordinates. The second is (u,r), where the retarded- 
time coordinate u is defined by u = t — r — 2M \n(r / 2M — 
1). The third is (v,r), where the advanced-time coordi- 
nate v is defined by v — t + r + 2M\n(r /2M — 1). In these 
coordinates the Schwarzschild metric takes the form 

g ab dx a dx b = f dt 2 + f- x dr 2 , (2.2) 
= fdu 2 -2dudr, (2.3) 
= fdv 2 + 2dvdr, (2.4) 

where / := l — 2M/r and M is the mass of the black hole. 
These systems share the property that the scalar r is 
adopted as one of the coordinates. Our formalism is not, 
however, limited to these coordinate choices; one retains 
the freedom of using any coordinate system whatever, for 
example, harmonic coordinates, isotropic coordinates, or 
double-null coordinates. 

We introduce the dual vector 



dr 



(2.5) 



which is normal to the surfaces of constant r(x a ); in the 
coordinates of Eqs. r a = (0, 1). We use g , 

the inverse to g a b, to raise its index: r a = g ab rb. This 
allows us to give a covariant definition to the function / 
that appears in Eqs. (|2.2|l - (|2.4() : 



2M 

/ = ! • 

r 



(2.6) 



We also introduce e a bi the (antisymmetric) Levi-Civita 
tensor on M 2 ] in the coordinates of Eqs. (|2.2II 12. 41) . 
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e tr = s ur — £vr = 1- The timelike Killing vector of the 
Schwarzschild spacetime is tangent to M 2 and is given 
by 



-e ah n; 



(2.7) 



in the coordinates of Eqs. t<X = (M)- We 

have t a i a = — / and i a r a = 0, and the vectors r a , t a 
form a basis on M 2 . In terms of this basis we have 

g ab = f-l(_ t a t b + r a r fe) and £ ab = _j-l^a r b _ r a f by 

The covariant derivative operator compatible with g ab 
is denoted V a ; we thus have V a gbc = 0. It is easy to 
show that for the Schwarzschild solution, 



V V & r 



M 

~~^9ab. 



(2.8) 



so that O = 2M/r 2 , where □ := g ab V a V b is the Lapla- 
cian operator on M 2 . We also have V a t b = (M/r 2 )e ab , 
which confirms that t a is a Killing vector. The Riemann 
tensor on M 2 is R abcd = (2M/r 3 )(g ac g bd - g a d9bc)- 

We let n AB be the inverse to Hab, the metric on the 
unit two-sphere. The covariant derivative operator com- 
patible with O ab is denoted Da; we thus have Da^bc = 
0. The Levi-Civita tensor on the unit two-sphere is de- 
noted sab, and sg^, = sin 9. The Riemann tensor on the 
unit-sphere is Rabcd = Qac^bd - VLad^bd- 

Covariant differentiation in the Schwarzschild space- 
time can be defined in terms of covariant differentiation 
in the submanifolds M 2 and S* 2 . If T a bc is the connection 
associated with V , and if T BC is the connection asso- 
ciated with Da, then it is easy to show that the nonvan- 
ishing components of the spacetime connection are given 

hv 4 F a — r a 4 F a — —rr a Oo^ i V A — r~ 1 r t) A 
oy 1 be — 1 he 1 BC " S 'BC, 1 Be ' 'c° B> 

and 4 T A BC = T A BC . Using these rules we find that the 
conservation identities for a stress-energy tensor T a @ in 
the Schwarzschild spacetime take the form 



W b T ab + D B T aB + -r b T ab ~ rr a n AB T AB = (2.9) 
r 



and 



V a T aA + D B T AB + -r a T aA = (2.10) 
r 



when expressed in terms of the submanifold connections. 



Vectorial spherical harmonics come in two types. The 
even-parity harmonics are defined by 



\rlra 
Y A 



D A Y ln \ 

while the odd-parity harmonics are 



X l r := -e?D B Y u 



(3.1) 



(3.2) 



Their components are listed explicitly in Appendix El 
The vectorial harmonics satisfy the orthogonality rela- 
tions 



and 



Y x tYA m ' dn = l(l + l)S w S n 



x A m x l A m ' dn = i(i + i)5 w s n 



(3.3) 



(3.4) 



in which an overbar indicates complex conjugation and 
dQ := sinOdOdcf) is an element of solid angle. We also 
have 



\rA yl'r. 
I lm^-A 



dn = o, 



(3.5) 



which states that the even-parity and odd-parity har- 
monics are always orthogonal. The definitions l|3.1|l and 
(|3.2|l for the vectorial spherical harmonics are identical 
to those provided by Regge and Wheeler 0. 

Tensorial spherical harmonics come in the same two 
types. The even-parity harmonics are SIabY 1 ™ 1 and 



Ylm 
AB ■ — 



D A D B + -1(1 + l)Sl AB 



Y l 



while the odd-parity harmonics are 



x ab '■— ~-^{ s aDb 



£r D a)Dc:Y 



(3.6) 



(3.7) 



Their components are listed explicitly in Appendix ^ 
The tensorial harmonics satisfy the orthogonality rela- 
tions 



y ab y a ™' dn 



1 



(l-l)l(l + l)(l + 2)S lv 6 mm> (3.8) 



and 



III. SPHERICAL HARMONICS 

In this section we introduce the scalar, vector, and 
tensor spherical harmonics that are used in the decom- 
position of the metric perturbation. The tensorial nature 
of the spherical harmonics refers to the unit two-sphere, 
and in this section we use the metric O^g and its inverse 
n AB to lower and raise all upper-case Latin indices. All 
tensorial operations (including covariant differentiation) 
shall refer to this metric. 

The scalar harmonics are the usual spherical-harmonic 
functions Y lm (9 A ). They satisfy the eigenvalue equation 
[Q AB D A D B + 1(1 + l)]Y lm = 0. 



X h^ X AB d^l 



We also have 



~(l-l)l(l + l)(l + 2)5 

ll'Umm' • 

(3.9) 



Y AB x l ' A i dn = o 



and 



CiAB \rlm 
iL 1 AB 



o = n AB x l XB. 



(3.10) 



(3.11) 



The definition (|3.6|) for the even-parity harmonics does 
not agree with that of Regge and Wheeler |l| , who work 
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instead with the set n AB Y lm and D A D B Y lm . We find it 
more convenient to form the tracefree combinations Y^g , 
which have the property of being (pointwise) orthogonal 
to n AB Y lm . The definition for the odd-parity har- 
monics also differs from Regge and Wheeler's, but only 
by an overall minus sign, which we find convenient to 
introduce. 

The tensorial harmonics Y^g and X l AB can be related 
to the spherical-harmonic functions of spin-weight s — 
±2 [l6|. and to the pure-spin harmonics used by Thorne 
These relations are explored in Appendix 1X1 

IV. EVEN-PARITY SECTOR 
A. Perturbation fields and gauge transformations 

The even-parity sector refers to those components of 
the metric perturbation that can be expanded in terms of 
the even-parity spherical harmonics Y lm , Y A m , fl AB Y Lm , 
and Yj^. Introducing the notation 4 g a i, = g a b + p a b, 
i g a B = PaB, and 4 g AB = r 2 Vt AB +p A B for the perturbed 
metric, the even-parity sector of the metric perturbation 
is 

Pab = (4.1) 
I in 

PaB = £jl m Y<™, (4.2) 
I in 

Pab = r^iK^ABY 1 " 1 + G lm Y$£). (4.3) 
i ni 

In this section the sums over / are taken to begin at 
1 = 2; the low multipoles (1 = and I = 1) will 
be considered separately in Sec. IVIIII The fields h l *£, 
K lm 7 and G lm are defined on M 2 and they de- 
pend on the coordinates x a only. They are closely re- 
lated to the quantities first introduced by Regge and 
Wheeler pj , who worked exclusively in terms of the usual 
Schwarzschild coordinates (t,r). In these coordinates 
(discarding for brevity the spherical-harmonic labels) we 
have h tt = fH , h tr = Hi, h rr = H 2 /f, jt = h , and 
j r = hi. The function G introduced in Eq. (|4.3[l is iden- 
tical to the corresponding Regge- Wheeler quantity, but 
K is different: inhere = ^rw — hKl + the difference 
originates from the fact that Regge and Wheeler work 
with Q. AB Y lm and D A D B Y lm instead of Q AB Y lm and 

1 AB- 

Even-parity gauge transformations are generated by a 
dual vector field S Q = (S a ,EU) that is expanded as 

lm 

Ea = ^ImyUn. (4.5) 
lm 

the fields ^ m and £} m depend on the coordinates x a only. 



Under such a transformation the perturbation quantities 
change according to (see Appendix IB|l 

Kb -> Kb = h ab - V a £b - V 6 ^ a , (4.6) 

3a ~> 3a = 3a - & - V Q C + ^r a f, (4.7) 

K K> = K+ l -V±Vz-*r a Z a , (4.8) 

G G' = G-^, (4.9) 

where we have discarded the spherical-harmonic labels 
for brevity (we shall continue with this practice until the 
end of the section). It is easy to show that the combina- 
tions 

Kb ■= h ab - V a e 6 - Vfc£ Q (4.10) 

and 

K:=K+]-l{l + l)G--r a e a (4.11) 
2 r 

are gauge invariant, where 

£ a :=ia-^r 2 V a G. (4.12) 

Equations Ij4.7|l and (|4.9|l reveal that one can always 
choose a gauge in which j a = = G; this is the Regge- 
Wheeler gauge. Equations (|4.10() - (|4.12|l imply that h a \, = 
h ab and K = K in the Regge- Wheeler gauge. 



B. Perturbation equations 

The Ricci tensor of the Schwarzschild spacetime van- 
ishes, and as a consequence its perturbation is gauge in- 
variant. Its computation can therefore be carried out in 
any convenient gauge, and the Regge- Wheeler gauge is 
clearly convenient. The steps involved are as follows. We 
substitute Eqs. 14. ip 14.3fl . having set j a = G = 0, into 
the Ricci tensor of Appendix El We simplify the result 
and find that SR ao is expanded in terms of Y lm , SR aB 
in terms of Y B m , and SR AB in terms of both £l AB Y lm 
and Y l A 71 B . Fr om the Ricci tensor we compute the Ein- 
stein tensor and set the result equal to 8nT a ^ , which 
describes the material source for the gravitational per- 
turbations. Each spherical-harmonic component of the 
field equations can then be extracted by involving the 
orthonormality relations l|3.3f) and 13. 8|) satisfied by the 
spherical harmonics. At the end of this calculation we 
take advantage of the fact that h a b = h a b and K = K 
in the Regge- Wheeler gauge. This allows us to make 
the substitutions h a b — ► hab, K — > K and therefore to 
express the gauge-invariant Einstein tensor in terms of 
gauge-invariant quantities. 

Our final results are 
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;<ib 



Qa 



V c V {a h c b) - ^g ab V c S7 d h cd - ~0(h ab - g ab h) - ^V a V b h + ^r c (V \ a h c b) - g ab V d h cd ) 

'1(1 + 1) 2M~ 



- ^r c V c (h ab - g ab h) + l -^^— h ab - -^g ab r c rdh cd - 

2 3 (I — 11(7 4- 21 

- X7 a V b K + g ab DK - -r {a V b) K + -g ab r c V C K - K - ^ >-g ab K, 

V c h c a - V a h + -r a h - V a K, 
r 

Oh - V a V b h ab - -r a V b h ab + -r a W a h - „ , ft + OK + -r a V a K, 



g ab h 



2r 2 



= -h, 



(4.13) 
(4.14) 

(4.15) 
(4.16) 



where h := g h ab and □ :— g V a V b . The source terms 
are 



Q a = 



8tt J T ab Y lm dn, 
167rr 2 



1(1 + 1) 



T aByhn ^ 



Q b = 8nr 2 



T AB n AB Y lm dQ, 



(4.17) 
(4.18) 
(4.19) 



327rr 4 



(l-l)l(l + l)(l + 2) 



T AB YX% dQ. (4.20) 



In Eqs. H4.17fl - H4.20fl the stress-energy tensor is imag- 
ined to be given in fully contravariant form; T ab , T aB , 
and T are then its components in the spacetime co- 
ordinates (x a ,9 A ). In the event where the stress-energy 
tensor would be given in covariant or mixed form, its 
indices would have to be raised with (g ab , r~ 2 Q AB ) - 
the inverse Schwarzschild metric — before evaluating the 
source terms. In Eqs. ()4.13(l - (|4.16|l all lower-case Latin 
indices are lowered and raised with g ab and g ab , respec- 
tively. 

The perturbation equations are not all independent. 
By virtue of the Bianchi identities, or the conservation 
equations (|2.9|) and (|2.1U|I . they are related by 



ab 



-r b Q 



ab 



Ht ' -1 ).Q« -l r a Q b =0 (4.21) 



2r 



and 



VaQ a 



2r 2 W 



Q b = 0. (4.22) 



C. Master equation 

The Zerilli-Moncrief function is defined by 



im 



2r 



1(1 + 1) _ 



2 



K lm + ±_( r a r b~ h lrn _ ^ay^imj 



(4.23) 



where 



A:= (l-l)(l + 2) + — . 

r 



(4.24) 



This is a covariant and gauge-invariant generalization of 
the definition provided by Lousto and Price [18j. who 
worked with the usual Schwarzschild coordinates and in 
the Regge- Wheeler gauge. Their even-parity function 
agrees, up to a normalization factor, with the gauge- 
invariant function first introduced by Moncrief [9j . Mon- 
crief's function, in turn, is a variant of Zcrilli's original 
even-parity function [j| (Zerilli worked in the frequency 
domain instead of the Schwarzschild time domain). Up to 
the same normalization factor the Moncrief and Fourier- 
transformed Zerilli functions are equal up to the pres- 
ence of source terms; they are equal only where there is 
no matter. The normalization adopted in Eq. (|4.23|) will 
be seen to be convenient when we discuss gravitational 
radiation at future null infinity fSec, IVl|) and at the hori- 
zon (Sec. IVIIf) : our definition of the even-parity master 
function is well adapted to the description of radiation 
fields. 

The perturbation equations (|4.13|) - (|4.16[1 give rise to 
a wave equation for the function 'E'even (we resume our 
practice of omitting the spherical- harmonic labels). The 
manipulations are long and tedious, and we shall not 
present them here (the reader may consult Ref. [l4| for 
details). We simply state the final result: As a con- 
sequence of the field equations, the even-parity master 
function satisfies the Zerilli equation 



(I— 1 l / cvcn)^'cvcn — S even: 



(4.25) 



with potential 



1 



A 2 

and source term 



H + 2 6M 



36M 2 / 2M\ 



(4.26) 



6 



q - 4 r n a 1 ni I 2 

'-'even — i 1 a^ ' r , ^\ a 

A r (// + 2) A 



-2rVV a ( 



24M 
A 



r a r b Q ab 



M M 2 * 
- 2) + I2(fi - 3)— + 84— 



(4.27) 



where /x := (I - l)(l + 2) and Q := g a6 Q afc . The validity 
of Eqs. (|4.25|) - (|4.27|) can be verified by brute-force eval- 
uation of both sides of Eqs. Q4.25p . The general source 
term for the Zerilli equation has never been presented in 
the literature; we display it here for the first time, in a 
form covariant under an arbitrary transformation of the 
coordinates x a used in the two-dimensional submanifold 
Ai 2 of the Schwarzschild spacetime. 

V. ODD-PARITY SECTOR 



where we have discarded the spherical-harmonic labels 
for brevity (we shall continue with this practice until the 
end of the section). It is easy to show that the combina- 
tions 

K =h a - \v a h 2 + -r a h 2 (5.7) 
2 r 

are gauge invariant. Equation (|5.6(l reveals that one can 
always choose a gauge in which h 2 — 0; this is the Regge- 
Wheeler gauge. Equation (|5.7|l implies that h a = h a in 
the Regge- Wheeler gauge. 



A. Perturbation fields and gauge transformations 

The odd-parity sector refers to those components of the 
metric perturbation that can be expanded in terms of the 
odd-parity spherical harmonics X 1 ™ and X^jb . Recalling 
the notation 4 g ab = g ab + p ab , 4 g aB = p aB , and i g AB = 
r 2 flAB + pab for the perturbed metric, the odd-parity 
sector of the metric perturbation is 



Pab 
PaB 

PAB 



= o, 



Ei 1 7ii vim 

Irn 



hTx 



AB- 



(5.1) 
(5.2) 

(5.3) 



In this section the sums over I are taken to begin at / = 2; 
there is no odd-parity perturbation with I = 0, and the 
case I — 1 will be considered separately in Sec. IVIIII 
The fields h l ™ and h\ m arc defined on M 2 and they 
depend on the coordinates x a only. They are closely 
related to the quantities first introduced by Regge and 
Wheeler , who worked exclusively in terms of the usual 
Schwarzschild coordinates (t,r). In these coordinates 
(discarding spherical-harmonic labels) we have h t = ho 
and h r = h%. Except for a sign — refer back to the dis- 
cussion following Eq. (|3.11|) — the function hi is identical 
to the corresponding Regge- Wheeler quantity. 

Odd-parity gauge transformations are generated by a 
dual vector field S Q = (0, H^) that is expanded as 



"A 



lm 



Im \rlm 



(5.4) 



in which £ im depends on the coordinates x a only. Under 
such a transformation the perturbation quantities change 
according to (see Appendix lB|l 



h a 

h 2 



h' a = h a - V Q f + -r a £, 
r 



hi, 



2?, 



(5.5) 
(5.6) 



B. Perturbation equations 

The Ricci tensor of the Schwarzschild spacetime van- 
ishes, and as a consequence its perturbation is gauge in- 
variant. Its computation can therefore be carried out in 
any convenient gauge, and as in the preceding section 
we shall adopt the Regge- Wheeler gauge. We substitute 
Eqs. H5.1(l - H5.3(l . having set h 2 — 0, into the Ricci tensor 
of Appendix IbI We simplify the result and find that SR ab 
vanishes, 8R a s is expanded in terms of X 1 ^ 1 ', and SRab 
in terms of X 1 ^ . From the Ricci tensor we compute the 
Einstein tensor and set the result equal to 8irT a P . Each 
spherical-harmonic component of the field equations can 
then be extracted by involving the orthonormality rela- 
tions l|3.4(l and l|3.9() satisfied by the odd-parity harmon- 
ics. At the end of this calculation we take advantage of 
the fact that h a = h a in the Regge- Wheeler gauge. This 
allows us to make the substitution h a — > h a and therefore 
to express the gauge- invariant Einstein tensor in terms of 
gauge-invariant quantities. 

Our final results are 



Pa = 
P = 

where □ := 

pa _ 

P = 



-nh a + V a V b h b + -(r b V a h b 
r v 



r a V b h b ) 



1(1- 



- ^2 r ar b h b 
Vah a , 

- 5 afc V a Vb. The source terms are 

167rr 2 
1(1 + 1) 



(5.8) 
(5.9) 

(5.10) 



167rr 



(l-l)l(l + l)(l + 2) 



T AB X%dtt. (5.11) 



In Eqs. (|5.10() and 15.11(1 the stress-energy tensor is imag- 
ined to be given in fully contravariant form; T aB and 
T AB are its relevant components in the spacetime coor- 
dinates (x a ,9 A ). In the event where the stress-energy 
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tensor would be given in covariant or mixed form, its 
indices would have to be raised with (g ab , r~ 2 Q AB ) - 
the inverse Schwarzschild metric — before evaluating the 
source terms. In Eqs. (|5.8|) and (|5.9|) all lower-case Latin 
indices are lowered and raised with g ao and g ab , respec- 
tively. 

The perturbation equations are not all independent. 
By virtue of the Bianchi identities, or the conservation 
equation p.lOfl . they are related by 



The validity of Eqs. (|5.14() - (|5.16|l can be verified by 
brute-force evaluation of both sides of Eqs. ((5.14(1 . The 
general source term for the Regge- Wheeler equation has 
never been presented in covariant form in the literature. 
We display it here for the first time, and note that in the 
usual Schwarzschild coordinates (and in the frequency 
domain), S^dd agrees with the source term presented in 
Eq. (18) of Jhingan and Tanaka 13]. 



v„ P . + ; f . P . _ ('-!)('+ 2> P=0 . 



(5.12) 



D. Regge- Wheeler function 



by 



C. Master equation 



The Cunningham- Price- Moncrief function is defined 



Ira 
odd 



2r 



(i-m + 2) 



(5.13) 



where £ a & is the Levi-Civita tensor on the submanifold 
M 2 . Apart from a different normalization factor, this 
is a covariant generalization of the definition provided 
by Cunningham, Price, and Moncrief 0] i who worked 
with the usual Schwarzschild coordinates. Our definition 
(and normalization) agrees with the odd-parity function 
considered by Jhingan and Tanaka 0] , who also worked 
with the Schwarzschild coordinates, but in the frequency 
domain. The normalization adopted in Eq. (|5.13|l will 
be seen to be convenient when we discuss gravitational 
radiation at future null infinity fSec. lVI|) and at the hori- 
zon (Sec. IVII(1 ; our definition of the odd-parity master 
function is well adapted to the description of radiation 
fields. It is noteworthy that it can also be expressed as 



2r 



odd 



(i-m- 



i ah (d a l 



2 , 
-r a h 
r 



in terms of the original (gauge-dependent) perturbation 
quantities, and in terms of the partial differentiation op- 
erator; the function is nonetheless gauge-invariant and a 
scalar. 

The perturbation equations 1(5.8(1 and (|5.9J) give rise to 
a wave equation for the function ^ dd (we resume our 
practice of discarding the spherical- harmonic labels) . As 
in the preceding section we simply state the final result: 
As a consequence of the field equations, the odd-parity 
master function satisfies the Regge- Wheeler equation 



(□ - Kdd)*odd = S Q 



dd j 



with potential 



V 



dd 



1(1 + 1) 6M 



and source term 

•Sodd 



2r 



(/-!)(/ + 2) 



e ab V a P b 



(5.14) 



(5.15) 



(5.16) 



The Cunningham-Price-Moncrief function is a close 
cousin to the more familiar Regge- Wheeler function Q] , 
whose covariant and gauge-invariant definition is 



RW 



-r a h l : 

r 



(5.17) 



As we shall see in Sees. IVll and lVllI the Regge- Wheeler 
function is not well suited to describe the gravitational 
radiation field, and in this paper we adopt the function 
^odd °f Eq. (|5.13(l as the fundamental odd-parity master 
function. It is straightforward to use the perturbation 
equations to show that these functions are related by 



1 



RW 



dd 



(i-m + 2) 



-r n P\ 



(5.18) 



where t a — —e ab rb is the Killing vector of Eq. 1)2.7(1 . 
Outside of sources, and apart from a factor of one-half, 
the Regge- Wheeler function is the time derivative of the 
Cunningham-Price-Moncrief function. 

The function of Eq. H5.17|l also satisfies the Regge- 
Wheeler equation, 

(□ - Kdd^RW = Srw, (5-19) 

with the potential of Eq. (|5.15|l and a new source term 
given by 



RW 



r a (V a P-P a ) 



1 



P 



(5.20) 



This covariant form for the source term is also a new 
result. 



VI. RADIATION AT FUTURE NULL INFINITY 

To examine the gravitational perturbations near future 
null infinity we adopt the retarded coordinates (u, r, 0, (f>) 
and express the two-dimensional Schwarzschild metric in 
the form of Eq. 12.3(1 . In these coordinates, future null 
infinity corresponds to taking the limit r — > oo keeping u 
fixed, and our strategy will be to expand the metric per- 
turbations in powers of r . In asymptotically Cartesian 
coordinates the radiative part of the metric would scale 
as r , transforming to spherical coordinates produces 



the scalings 



T) rad 

Pab 



0{r~ 



T) lad 

PaB 



0(r°), 



Ti rad 

Pab 



0(r) (6.1) 
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for the radiative part of the metric perturbations. Our 
goal is to isolate this, and to calculate how much energy 
and angular momentum is carried away by the radiation. 
We find it advantageous to work in gauge in which 



t a p ab = = t a Pal 



(6.2) 



where t a is the Killing vector of Eq. (|2.7|) . In spite of 
the fact that t a is not a null vector (except on the even 
horizon, which is well outside our domain of considera- 
tion), this gauge happens to be a perfectly respectable 
"radiation gauge." 

In this section we take I > 2; as is well known, the 
low multipoles / = and Z = 1 do not contain radiative 
degrees of freedom. We assume that the matter distribu- 
tion responsible for the radiation is confined to a bounded 
volume, and that our domain of consideration is outside 
this volume; we shall therefore be solving the vacuum 
field equations. 

We begin with the even-parity sector of Sec. IIVI The 
gauge conditions imply h uu = h ur = j u = 0, and the 
scalings of Eq. (|6.1|) imply that we are looking for the r _1 
part of h rr , K, and G, as well as the r° part of j r . These 
can be determined by expanding the components of the 
metric perturbation in powers of r _1 and substituting 
them into the field equations of Eqs. H4.13j) - I|4.16[) . The 
expansions for h rr , K, and G begin at order r" 1 and 
each coefficient is a to-bc-dctcrmincd function of u; the 
expansion for j r begins instead at order r°. The field 
equations return Q a b, Qa, , and Q" expanded in powers 
of r _1 , and setting each coefficient to zero determines the 
metric perturbation. We obtain 



1(1+1) 



a(u')du' + 0{r- 4 ), (6.3) 



Jr 

K 
G 



a(u) 



r 

1(1 + 1) 
2r 3 



+ 0(r- 2 ), 



(6.4) 



a(u')du' + 0(r- A ), (6.5) 



+ 2) r 



0(r 



(6.6) 



where a(u) is a function that is not determined by the 
vacuum field equations, and a := da/ du. We see that the 
radiative part of the perturbation is contained entirely 
in the function G(u,r). With the metric perturbation 
of Eqs. H6.3(l - (|6.6|l we may evaluate the Zerilli-Moncrief 
function of Eq. H4.23[) . The result is 



(/-!)(/ + 2) 



a(u) 



Oir- 1 ), 



(6.7) 



and we conclude that the radiative part of the even-parity 
sector is given by 



We turn next to the odd-parity sector of Sec. El The 
gauge conditions imply h u — 0, and the scalings of 
Eq. 1)6. Q imply that we are looking for the r° part of 
h r and the r 1 part of These are determined by fol- 
lowing the same procedure as in the even-parity case, and 
we obtain 



h r 
hi 



(*-!)(* + 2) 
2r 

b(u)r + 0(r°), 



b(u')du' + 0(r~ 2 ), (6.9) 
(6.10) 



where b(u) is not determined by the vacuum field equa- 
tions. We see that the radiative part of the perturbation 
is contained entirely in the function ti2(u, r). With the 
metric perturbation of Eqs. (|6.9|l . i|6.10[) we may evaluate 
the Cunningham-Price-Moncrief function of Eq. H5.13|l . 
The result is 



* odd = b(u) + 0(r- 1 ), 



(6.11) 



and we conclude that the radiative part of the odd-parity 
sector is given by 



rad. 

Pab 



>dd 



■£*Sdd(«,r = oo)A: 



lm 
AB- 



(6.12) 



It is obtained by integrating the Regge- Wheeler equation 
(15.141) and evaluating the Cunningham-Price-Moncrief 
function at r = oo. Notice that by virtue of Eq. H5.18J) . 
the radiative field could instead be expressed in terms of 
the u-integral of the Regge- Wheeler function. The need 
to perform this integration is inconvenient, and it is the 
simple relationship of Eq. H6.12fl that has motivated the 
adoption of the Cunningham-Price-Moncrief function as 
the fundamental odd-parity master function. 

The full radiative field is obtained from Eqs. (|6.8|l and 
We have pl% d = = p^g and 



rad 

Pab 



l m 



lm 



rim 



*cvcn J AB 



ijrlm vim \ 
W odd A AS ) : 



(6.13) 



\\,lm . 
^cvcnV"' ' 



oo) and ^ 



lm 
dd 



where ^ 

oo). As expected, the radiative field is transverse, and 
tracefree by virtue of Eqs. (|3.11(l . The two fundamen- 
tal polarizations of the gravitational wave can be defined 



by h+ := pee/r 2 and h > 



P8rp/(r 2 sin 9). Using the 



components of the tensorial spherical harmonics listed in 
Appendix^ we obtain 



lm 
even 



odd" 



d0 2 

d_ 

86 ' 



+ -1(1 + 1) 



Y 



im 



COSf 



(6.14) 



and 



rad, even \ ^ , T ,; 

Pab = r l^ *ov 01 , 

lm 



(u, r = oo)Y j 



l m 
AB- 



(6.8) 



It is obtained by integrating the Zerilli equation (|4.25|) 
and evaluating the Zerilli-Moncrief function at r = oo. 



r / y 1 even. 



lh 



sm ( 



00 



sinf 



Y 



! iii 



odd 



1) 



(6.15) 
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The energy and angular momentum carried away by 
the gravitational radiation can be calculated using the 
techniques developed by Thorne |l7j |. We compare our 
Eq. (|6.13l) to his Eq. (4.3), taking into account the rela- 
tionship between our tensorial harmonics and his pure- 
spin harmonics (this is spelled out in Appendix 
Thome's mass multipole moments are thus seen to be 
proportional to *if™ en (u, r = oo), and his current mo- 
ments are proportional to ^^SAu,r = oo). Substituting 
these into Thome's Eq. (4.16) we obtain 



64tt 



lm 



odd 



(6.16) 



for the averaged rate at which the energy escapes to 
future null infinity. Substituting instead into Thome's 
Eq. (4.23) returns 



dJ_ 

du 



1 

64^ 



^2(1- 1)1(1 + l)(l + 2){im) 



I in 



\U L7Tl \T/ im 

even even 



^ odd odd 



(6.17) 



We assume that the matter distribution responsible for 
the radiation does not come near the event horizon; we 
shall therefore be solving the vacuum field equations in 
an empty neighborhood of the event horizon. 

We begin with the even-parity sector of Sec. IIVI The 
gauge conditions of Eq. I|7.1|) imply that h vv = h vr = 
j v = 0, so that h rr , j r , K, and G are the only non- 
vanishing components of the metric perturbation. As 
stated above, the radiation crossing the event horizon is 
described entirely by K and G evaluated at r = 2M. The 
field equations, however, imply that K vanishes on the 
horizon. This is verified by expanding h rr , j r , K, and G 
in powers of / = 1 — 2M/r and substituting them into 
the field equations of Eqs. (|4.13|l " (|4.16|l . This calcula- 
tion reveals also that G(v,r = 2M) = a(v), a function 
that is not determined by the vacuum field equations. 
On the other hand, the expansions allow us to evaluate 
the Zerilli-Moncrief function of Eq. I|4.23|l . and the result 
is \l/ even (r = 2M) — 2Ma(v). We conclude that on the 
horizon, the radiative part of the even-parity sector is 
given by 



lm 
even 



{v,r = 2M)Y A m B - (7.2) 



for the averaged rate at which the angular momentum 
escapes to infinity. This is the component of the angular- 
momentum vector in the arbitrary z direction which de- 
fines the orientation of the angles 8 and cf>. The overbar 
indicates complex conjugation, and it is not difficult to 
show that (dJ/du) is real. The averaging carried out in 
Eqs. (|6.16|l and (|6.17|) is over a characteristic time scale 
associated with the gravitational wave. 



VII. RADIATION AT THE EVENT HORIZON 

To examine the gravitational perturbations near 
the event horizon we adopt the advanced coordinates 
(v, r, 0, 4>) and express the two-dimensional Schwarzschild 
metric in the form of Eq. I|2.4I) . We want to calculate how 
much energy and angular momentum is transfered to the 
black hole by the perturbation, and we shall do so by 
following the methods devised by Poisson 0] . As in the 
preceding section we impose the gauge conditions 



t a p ab = = t a PaB 



(7.1) 



on the metric perturbation; t a is still the Killing vector 
of Eq. (|2.7(l . Recall that this vector is null on the event 
horizon, and Eq. I|7.1|l . evaluated at r = 2M, is equiva- 
lent to Poisson's Eq. (6.5). Poisson then shows that the 
part of the metric perturbation which is associated with 
the transport of energy and angular momentum across 
the horizon is pab- [Refer to Poisson's Eq. (6.10), which 
establishes the equality between pab and ^\ B ^ the per- 
turbation of the horizon's intrinsic metric] 

In this section we take I > 2, as the low multipoles I = 
and I = 1 do not contain radiative degrees of freedom. 



It is obtained by integrating the Zerilli equation (|4.25() 
and evaluating the Zerilli-Moncrief function at r = 2M. 

We turn next to the odd-parity sector of Sec. El The 
gauge conditions imply h u — 0, so that h r and /12 are 
the only nonvanishing components of the metric pertur- 
bation. The radiation crossing the event horizon is de- 
scribed entirely by /12 evaluated at r = 2M, and the 
field equations imply that b(v) := h,2(v,r — 2M) re- 
mains as an undetermined function. They also imply 
that h r (v,r = 2M) = c(v), with 



c(v) 



8M 2 



b(v') dv' . 



These statements are verified by expanding h r and Y12 
in powers of / = 1 — 2M/r and substituting them 
into Eqs. I|5.8|l and (|5.9|l . The expansions allow us also 
to evaluate the Cunningham-Price-Moncricf function of 
Eq. H5.13fl . and the result is 



* odd (r = 2M) 



4M 



dc 



(l-l)(l + 2)dv 



1 

2M 



b(v). 



We conclude that on the horizon, the radiative part of 
the odd-parity sector is given by 

pTb ° dd = -2M J2 * odd («, r = 2M)X% . (7.3) 

lm 

It is obtained by integrating the Regge- Wheeler equation 
(|5.14l) and evaluating the Cunningham-Price-Moncrief 
function at r = 2M. Notice that by virtue of Eq. I|5.18|l . 
the radiative field could instead be expressed in terms of 
the w-integral of the Regge- Wheeler function; this would 
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give rise to Poisson's Eq. (7.3), a less convenient expres- 
sion. 

The full radiative field is obtained by adding Eqs. Q7.2D 
and J72J|. We have 

p% = 2M (*M3 V^X'Tb) , (7-4) 

lm 

where ¥™ en = ¥™ en (v, r = 2M) and = ¥% d (v, r = 
2M). This result should be compared with Eq. H6.13JI . 
The rates at which the gravitational perturbation trans- 
fers energy and angular momentum to the black hole can 
now be calculated using the method described in Sec. VII 
of Poisson pjj. Our Eq. (|7.4|) replaces his Eq. (7.5), and 
the rest of the calculation is identical. The final results 
are 

(§) = i£«-w< + w + 2> 

v lm 

x(|^ v "c„f +|^X i f) (7.5) 

and 

(%) = ebp-iW + i)(i + 2)«m) 

x ' lm 

x (¥ m i> lm +¥ 7 !l,& n } A ). (7.6) 

\ ^ even even ' odd odd / V 1 w / 

These equations replace Poisson's Eqs. (7.8) and (7.9). 
Notice the similarity between Eqs. I|7. 5fl and (|6.16l) . and 
between Eqs. <|7.6[) and l|6.17|) . In Eq. I|7.6|) . J represents 
the component of the hole's angular-momentum vector 
in the arbitrary z direction which defines the orientation 
of the angles 9 and </>. The overbar indicates complex 
conjugation. The averaging carried out in Eqs. I|7.5|l and 
(|7.6(l is over a characteristic time scale associated with 
the gravitational perturbation. 

VIII. LOW MULTIPOLES 

To conclude our presentation of the perturbation for- 
malism we now handle the special cases I = and 1 = 1. 
Our discussion here will essentially reproduce Zerilli's 
Appendix G 0, but we will frame it in the language 
developed in this paper. For concreteness, and for sim- 
plicity, we will focus mostly on the task of integrating the 
vacuum field equations for these low multipoles. 

A. Monopole perturbation 

When I = the only relevant spherical harmonic is 
V 00 , which is a constant. It follows that Ya = Xa = 
Yab = Xab = 0, and the perturbation is of even parity. 
The only nonvanishing metric perturbations are 

Pab = K b Y 00 , pab = r 2 Kn AB Y w ; (8.1) 



the fields h a and G are not defined. The freedom to 
perform a gauge transformation is contained in S a = 
£, a Y 00 , Sa = 0, and the perturbations transform as 

Kb -> h' ab = h ab - V a £& - V & £ a , (8.2) 

K -> K' = K- -r a £ a . (8.3) 

r 

There is no analogue here of the gauge-invariant quanti- 
ties h a i, and K that were introduced in the general case. 
The relevant field equations for h ab and K are the Q ab 
and Q b equations of Eqs. (|4.13() and (|4.15() . respectively, 
in which we set I = 0, h ab = h a b, and K = K. The Q a 
and equations of Eqs. (|4.14() and (|4.16() . respectively, 
are not defined. 

We now solve the perturbation equations in the ab- 
sence of material sources; in this case Q a b = = Q b . 
We adopt the advanced coordinates (v,r,8,<j>) and the 
metric of Eq. (|2.4|l . For convenience we rescale the value 
of F 00 and set it equal to unity. We choose £ a so as to 
enforce the gauge conditions h rr = K = 0. This does not 
fully determine the gauge, because these conditions are 
preserved under an additional transformation generated 
by £,v = —fa(v) and £ r = a(v), where / = 1 — 2M/r 
and a(v) is an arbitrary function. This transformation 
corresponds to a redefinition of the advanced-time coor- 
dinate v, of the form v — > v ' = v + = v + a{v), and 
it changes the nonvanishing perturbations according to 
h vv — > Kv = Kv + 2/a(u) and h vr — ► h' vr = h vr - a(v), 
where an overdot indicates differentiation with respect to 
v. 

We now involve the field equations. We notice first 
that the Q rr = equation implies that h vr is in fact in- 
dependent of r; since this function of v can be altered 
at will by a gauge transformation, we may set h vr = 0. 
The gauge is now fully determined, and the sole remain- 
ing perturbation field is h vv . By virtue of the Q vr = 
equation we have that rh vv must be a function of v only, 
but the Q vv = equation constrains this function to be a 
constant. We therefore have h vv = ( const ant )/r, and this 
completes the integration of the perturbation equations. 

The vacuum, monopole perturbation of a 
Schwarzschild black hole is given by 

Pw = , 8.4) 

r 

where SM is a constant; all other components vanish. 
This addition to g vv simply shifts the Schwarzschild mass 
parameter from M to M + SM. The perturbation must 
respect Birkhoff's theorem, and the perturbed metric is 
indeed another Schwarzschild solution. 



B. Odd-parity dipole perturbation 

The only surviving odd-parity spherical harmonics for 
I = 1 are X A m , which are obtained from V 1 ™ 1 using 
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Eq. (|3.2|l . The tensorial harmonics X^ vanish, and the 
only surviving components of the metric perturbation are 

PaB =Y, h a nX B m - (8-5) 

m 

The perturbations h^ m can be altered by a gauge trans- 
formation generated by S a = 0, Ea = J2 m £, lm X A m ; they 
change according to 

K^h' a = h a - V a £ + -r„£. (8.6) 
r 

There is no analogue here of the gauge- invariant fields h a 
introduced in Sec. El for the general case. The relevant 
field equation for h a is the P a equation of Eq. I|5.8|l . in 
which we set I = 1 and h a = h a . The P equation of 
Eq. I|5.9[> is not defined. 

We now solve the perturbation equations in the ab- 
sence of material sources; in this case P a = 0. We once 
more adopt the advanced coordinates (v,r,0,<p) and the 
metric of Eq. I|2.4|) . For convenience we replace the set 
of complex harmonics Y lm with the set Y m given by 

(Y°,Y S ,Y C ) = (cos 9, sin 9 sin <f>, sin 9 cos <j>) 

= : (ft 3 , ft 2 , ft 1 ). (8.7) 

This set is real, but not normalized; summation over m 
will be replaced by a summation over the abstract index 
m = (0, s,c). Notice that ft 1 represents the Cartesian 
components of a unit radial vector; in ordinary vectorial 
notation, £1 — r/r, where r = (x,y,z) and r = \r\. We 
shall also need ft^ := dQ t /d9 A , the derivatives of ft 1 with 
respect to the angles 9 A . 

We choose £ so as to enforce the gauge condition 
K = 0. This does not fully determine the gauge, because 
this condition is preserved under an additional transfor- 
mation generated by £ = a(v)r 2 , where a(v) is an arbi- 
trary function. This transformation changes h v accord- 
ing to h v — > h' v = h v — a(v)r 2 , where an overdot indicates 
differentiation with respect to v. 

We now involve the field equations. We note first that 
the P r = equation implies that h v must have the form 
of h v — ki(v)r 2 + &2 («)/?", where k\ and &2 are arbitrary 
functions. The first term can be removed by a gauge 
transformation, and we are left with h v = hi (v) jr and a 
fully determined gauge. The P v = equation then im- 
plies that ki must in fact be a constant. We therefore 
have h v = (constant)/r, and this completes the integra- 
tion of the perturbation equations. 

The vacuum, odd-parity, dipole perturbation of a 
Schwarzschild black hole is given by 

m 

where (SJ) m are the constants identified in the preced- 
ing paragraph, and X™ — —£ A B DBY m are the real odd- 
parity harmonics; all other components of the metric per- 
turbation vanish. Relating the spherical harmonics to 



the unit radial vector ft 1 as in Eq. (|8.7|l . and writing the 
summation over m as an implicit summation over i, we 
obtain our final expression, 

PvA = --e/(W)ift* B . (8.8) 
r 

To give an interpretation to this result we perform a 
transformation from the spherical coordinates (r, 6 ) to 
the Cartesian coordinates x 1 = rft*. The metric per- 
turbation transforms as p V i = r~ x p v A$l A , where Slf := 
n AB 5 %3 tt j B . This gives p m = -2r- 2 e AB tt A nf(SJ) j . We 
now involve the easily-established identity EAB^f^f = 
Eyfeft , where is the permutation symbol (a fully an- 
tisymmetric tensor with entries —1, 0, or 1). This gives 

Pvi = -^e ijk (6J) j fl k = ™(«J x r)i. (8.9) 

This form for the metric tensor allows us to identify 
8 J with the angular-momentum vector of the spacetime. 
The perturbed metric therefore describes a slowly rotat- 
ing black hole; it is the Kerr metric linearized with re- 
spect to its angular-momentum parameter. 



C. Even-parity dipole perturbation 

The even-parity spherical harmonics for I — 1 are Y lm , 
F4 1 ™, and VIabY 1 " 1 ] Yab vanishes. The metric pertur- 
bation is then 

Pab = J2 h lTY lm , (8.10) 

1 n 

PaB = (8.11) 

m 

Pab = r 2 ft AB ^^ lm r lm , (8.12) 

m 

and the fields G lm are not defined. Gauge transforma- 
tions are generated by S a = ^2 m £,l m Y lm and Ea = 
Em^ 1 " 1 ^]" 1 ' The perturbations change according to 

Kb -» h' ab = h ab - V a £& - VfeCa, (8.13) 
3a -» Ja = 3a -ia- V a £ + ^r a £, (8.14) 
K -» K' = K + - VCa- (8.15) 

There is no analogue here of the gauge-invariant fields h a b 
and K introduced in Sec. \N\ for the general case. The 
relevant field equations for h a b, j a , and K are the Q a b, 
Qa, Q b equations of Eqs. H4.13[) - l|4.15|> . m which we set 
I = 1, Kb = Kb - V a j b - Vbj a, and K = K — 2r~ 1 r a j a . 
The Q" equation of Eq. 14.16JI is not defined for I = 1 . 

We now solve the vacuum field equations, Q ao — Q a = 
Q b = 0. We adopt the advanced coordinates (v,r,9,(f>), 
the real spherical harmonics of Eq. I|8.7|l , and a gauge in 
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K = 0. This gauge is preserved under 



which h rr = / 
a transformation generated by 

r 

= ai(u), 
£ = a x (v)r + a 2 (v)r 2 , 

where a\{v) and az(v) are arbitrary functions. This has 
the following effect on the remaining components of the 
metric perturbation: 



h v 

Il r 



_ i , „. 2M . . , 2M 

Kv = h w - 2a 2 r — (2qi - a 2 ) + — 2~ a i; 



I r 



Jv 



Jv=Jv~ 



(di + a 2 ), 

.,2 



2M 

(«! + a 2 )r a\. 

r 



The Q rr = equation implies that h vr is a function of v 
only, and this component can be set equal to zero by a 
gauge transformation. The remaining gauge freedom is 
now restricted by a 2 — —a±. The Q r = equation im- 
plies that j v must be of the form j v = k\(v)r 2 + k 2 (v)/r, 
where k\ and k 2 are arbitrary functions of v. With this 
information and the Q v = equation we find that h vv 
must of the form h vv = — k 2 /r 2 + 3k 2 /r + ks + 2kir, where 
&3 is an additional arbitrary function. The Q vr = 
equation now reveals that fc 3 = 0, the Q vv = equa- 
tion yields 2Mk\ = —k 2 , and the task of integrat- 
ing the perturbation equations is completed. Defining 
k(v) := k 2 (v)/ (2M), we have obtained 

... , 6M ; . . 2M, . , 
h vv = -2k(v)r + ——k(v) —k{v) 



and 



w ^ 2 2M 
j v = —k(v)r H k{v) 



The vacuum, even-parity, dipole perturbation of a 
Schwarzschild black hole can be presented in a gauge in 
which the nonvanishing components are 

6M 2M 
p vv = -2r& t n l + djO 1 T aifl l (8.16) 



and 



PvA = -r 2 d. i ^ + ^^aiQ l A , 



2M 

T 



.17) 



where the vector a(v) is an arbitrary function of the 
advanced-time coordinate v. This perturbation, however, 
can be removed by a gauge transformation generated by 



r 

E A = ~r 2 aiVl l A + ra z n A . 



(8.18) 

(8.19) 
(8.20) 



The perturbed metric is therefore a Schwarzschild solu- 
tion expressed in a coordinate system which differs from 
the original system (v, r, ). These coordinates are not 
inertial. In the limit M the perturbed metric is 
recognized as the metric of flat spacetime expressed in 
a light-cone coordinate system centered on an acceler- 
ated world line (see Ref. [2(|); the acceleration vector 
is &i(v). For M ^ the perturbed metric describes a 
Schwarzschild black hole moving on "the same" acceler- 
ated world line. 
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for the nonvanishing perturbation fields. The remaining 
gauge freedom is contained in 

2M . . 

£„ = -a[v)r-\ atv), 

r 

£ r = a(v), 
£ = — a(v)r 2 + a(v)r, 

where we have set a := ot\ and a 2 = —a. Under this 
transformation the perturbations change according to 

, / , „ .. . \ 6M . , , 2M 
h m — ► vv = h vv + 2a{v)r — a(v) + —^-a(v) 



and 



Jv -> J v = Jv + a{v)r a{v). 



Setting a(v) — k(v) produces h' vv — j' v = 0, and we 
conclude that the perturbation is pure gauge. 



APPENDIX A: COMPONENTS OF THE 
VECTOR AND TENSOR SPHERICAL 
HARMONICS, AND RELATIONSHIP WITH 
OTHER SPHERICAL HARMONICS 

We first list the components of the tensorial harmonics 
introduced in Sec. IIIII According to Eq. (|3.1|l we have 



\rlm 



d_ 

00 

d 



-yira 



-Y 



According to Eq. 1)3. 2f) we have 

ylm _ 1 d 



sin 8 i 

d 



Y 



I III 



S m0^-Y lm . 
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According to Eq. 1|3.6[) we have 



im 



Y 



yim 
1 9<t > 



d 2 1 



Y 



! 111 



cos 9 d 



smf 



Y 



Irn 



Y 



hn 



d 1 

sin 9 cos 9— + -1(1 + 1) sin 2 i 



_d(j) 2 ' QQ ' 2 

And according to Eq. I|3.7() we have 



Y 



I in 



i 



■ (- 



d 2 



cos 9 d 



60 ~ sme\d9d(j) sin0. 



X$ = ^(sin 



Ai™ = sin 



i:) 2 



1 d 2 



d9 2 sin 9 d<f, 

u ~ -cosily'" 



89 J 



The tensorial harmonics F^g and X 1 ^ can be related 
to the spherical-harmonic functions of spin-weight s = 
±2 [l6| . Let m A and be a complex orthonormal 
basis on the unit two-sphere, with m A = 2~ 1 / 2 (1, i sin^). 
The relationship is then 



1 AB — 



- v /(Z-l)/(/ + l)(Z + 2) 
x (_ 2 Y lm m A m B + 2 Y lm m A rh B ) 



and 



x% = --v(i-i)i(i+m+ , 2) 

x (_ 2 Y lm m A m B - 2 Y lm m A fh B ), 

where s Y lm are the spin-weighted spherical harmonics. 
These equations can be compared with Eqs. (2.38e) and 
(2.38f) of Ref. [13. This reveals that our tensorial har- 
monics are intimately related to the "pure-spin" harmon- 
ics used by Thorne. The relationship is 



\rlm 
1 AB 



and 



^(i -1)1(1 + • 

^1(1- 1)1(1 +l)(l+ 



2)T 



E2,lr, 
AB 



2)T 



B2dm 



AB 



Notice that the pure-spin harmonics T^' lm and T^' m 
are normalized on the unit two-sphere. Our convention 
here differs from Thome's, who inserts a factor of in 
m A and fh A in order to normalize them on a two-sphere 
of radius r. 



9aB 

i 

9AB 



PaB, 



r 2 n 



AB + PAB , 



where (g a b, r 2 £l AB ) are the components of the back- 
ground metric and (PabiPaB,PAB) are the components 
of the perturbation. In this Appendix we allow the back- 
ground metric to be completely general, so long as it is 
spherically symmetric; we do not restrict it to be the 
Schwarzschild metric. We raise lower-case Latin indices 
with g ab , the inverse to g ab , and we raise upper-case Latin 
indices with Q AB , the inverse to AB . The inverse per- 
turbed metric is thus 



V 6 = g ab -p ab , 

4 aB 1 aB 

g = — 2P , 



4 AB 1 r,AB 1 AB 

9 = -7" - -TP , 



up to terms quadratic in the perturbations. 

Covariant differentiation with respect to the coordi- 
nates x a on the submanifold Ai 2 is indicated with V a : 
Vaffbc = 0. Covariant differentiation with respect to the 
coordinates 9 A on the unit two-sphere is indicated with 
D A : D A Vl B c = 0. Quantities which depend only on x a 
are covariantly constant relative to the connection T A BC ; 
for example D A r = 0. Quantities which depend only 
on 9 A are covariantly constant relative to the connection 
r a hc ; for example V a Q AB = 0. 

A straightforward calculation returns the components 
of the perturbed connection, which we denote + 
5T a 01 . We obtain 

Ar a — n a , 

01 be — ° fcc; 

sr a bC = \(d cP \ + v bP a c - v a PbC ) - ln P a c , 

ST a BC = ^(D BP a c + Dcp a B -V a PB c)+rr m n BC p am , 



ST A bc = ^(V b p A + V c p A ~D A p bc ), 



.•< - ^( D cPb A ~ D A Pbc + VbP A c ) ~ ^r bP A c , 



ST A BC = ^C A BC + l -r m n BC p mA , 



where r a := V a r, 



C a bc := ^cP a b + V b p a c -V a Pbc ), 



APPENDIX B: PERTURBATION OF THE RICCI 
TENSOR FOR A GENERAL 
SPHERICALLY-SYMMETRIC BACKGROUND 
SPACETIME 



and 



C A bc ■= ^(DcP A B + D BP A C ~ D A PBC ) 



The perturbed spacetime metric is written as 

gab = g a b+Pab, 



The Ricci tensor for the perturbed spacetime is equal 
to the background Ricci tensor R a p plus its perturbation 
8R a (3- We obtain 
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SRn.h — 



W m C m ab + ~r m C m ab ~ ^aV b p m m - ±D M D MPab + ±D M (V aPb M + W bPl! 
^V a V b p A { f + ^ (r a V bP M M + r b V a p M M ) - ^(r a r b - rV a V b r)p M M , 



2r 3 



1 



SRaB = ^D B [V m p m a - V a p m m + -r aP m m ) - -(Qp aB - V m V a p m B ) - -{r a V m p m B - r m y a p m B ) 



2 



-^(r a r m + rV a V m r)p m B + ^D M (D BPaM - D M p aB ) + ^^a{D M p M B ~ D B p M M ) 



r 
1 

2^ 



^r a (D M p M B -D B p M M ), 



SRab — &>AB 



1 



rr m V 



n(p mn - \g mn p\) + (r m r n + rV m V„r)p ron 



l -D A D BP m m + iv m (ZWB + D BP m A ) 



-r m £l AB D M p' 



i.M 



\ppAB + ^D M C M AB - —D A D B p a M 



1 



1 



-r m V m (p AB - -tt AB p 



M 
M 



2/1 



M 

M I J 



where □ := g ab V a V b . These expressions can be simplified by involving Eqs. (|2.6I) and (|2.8|l when the background 
spacetime is the Schwarzschild spacetime. 

Under a gauge transformation generated by the dual vector field S Q = (S a ,EA), the components of the metric 
perturbation change according to 

Pab — * Pab '■= Pab - V Q S & - V&S Q , 

, 2 

PaB — > P a B := PaB - V a ^, B - L>B^ a + ~r a ^ B , 

r 

Pab -* Pab : = Pas - -Da^b - D B E A - 2rr m E m tt AB . 



It can be shown that when the background Ricci tensor vanishes, 5R abl SR aB , and 5R AB are all invariant under this 
transformation. We use this property in Sees. IIVI and Ivl 



APPENDIX C: PERTURBATION EQUATIONS IN (t,r) COORDINATES 

In this Appendix the two-dimensional Schwarzschild metric is written in its standard form 

ds 2 = - / dt 2 + f- 1 dr 2 , 
with / = 1 — 2M jr. In these coordinates the wave operator is given by 

| 1 y fid 2 d ^ , _ 

D *=[-jdr 2 + d-r f d- r r- 

We set A := /(/ + 1) = fi + 2 and fx := (I - + 2) = A - 2. 

The even-parity perturbations are h tt , h tr , h rr , j t , j r , K, and G. Under a gauge transformation generated by £ t , 
£ r , and £ they change according to h u — > h' u — htt + A/i tt , etc., with 

Ah tt = -2— 6 + 

at r z 

9 , 9 , 2M _ 



<9r^ 1 9t^' ' r 2 / 
9 2Af 

d 



ar r 



AK = -y£r + ^£, 

AG = ~t 
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The gauge-invariant fields are 

htt 

htr 



2 ai 3 * 



£\2 

2-lr+r 2 —G-Mf — G, 
r z dt z or 



IMS 



d_ 

or" 



d_ 

di J 



■m 2 d 2 Q r-3M d c 



dt 



dr^ r r 2 f 



2M 

—Jr + r 



2 d 2 Q 2r-3M c 
dr 2 f dr 



K = 



r dr 2 



The field equations are 



Q tr 

Q* 
Q r 



3r-5M d - / <9 r (A + 2V + 4M- /x ~ 



dr 2 r 2 f dr 



d 2 



dtdr 



K + 



-W 2K + 



r-3M d T> f_d_~ ^\_~ 
~r 2 f~di ~ rdi rr ~ 2f tr ' 

(r-M)f d - 2/ d ~ fdr, Ar + 4M - / 2 , 



a " X + ~7^ htr TT htt + 

ar r at r Or 



2r 3 



^ti Trh rr — 



2r 2 



K. 



d ~ d - 19 - 2M- 

Tr h rr — — htr T" ~p "777-"- ^IHr, 

dt or j dt r z j 



ot ar ar 



r — M ~ (r - M)/ ~ 



ot^' 



9 2 - a 2 - 
otor dr z 



1 92 k 92 K 

f dt 2 dr 2 



2{r- M) d ~ r~3Md~ 
-htr T7 — -rrhtt 



r 2 f dt 



(r-M)fd- 2(r-M)d~ Ar 2 - 2(2 + X)Mr + AM 2 ~ 

-h rv H 7. ~7T~ K H ~ 1 ro htt 



r 2 f dr 
Ar 2 - 2/xMr - AM 2 - 



dr 



dr 



2r 4 / 2 



2r 4 



-J h u — fhrr- 



The conservation (Bianchi) identities read 





a 

at' 
a 



°=dt Qt 



a_ 

dr 
d 



Mf^„ 2r-hM 



= i-Q tr + ^-Q rr + —fQ tt + - 9 r" Q rr 



9r 



r 2 / 

r 2r^ 



2r 2 ^ r^' 



The Zerilli-Moncrief function is 



^even — 



2r 
+ 2 



A ' ^(./ 



■4* 

ar 



where A = /x + 6M/r, and its source term is 



S even - ^ ^ + (/i + 2 ) A 



^/|:(/o tt -7<r 



M M 2 
- 2) + 12(/i - 3)— + 84— 



24M 



Q rr + 2rfQ b 



fQ tt - - £ Q r 
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The odd-parity perturbations are ht, h r , and h%. Under a gauge transformation generated by £ they change 
according to h t — > h' t = h t + Ah t , etc., with 

Ah r = 

or r 
Ah 2 = -2f. 



The gauge-invariant fields are 



The field equations are 



ht = h t -- w h 2 , 

l t 13. 1, 

K = K - -TT«2 + -hi. 

2 or r 



2 ~ d 2 ~ 2d ~ \r-AM~ 
" — "a^a + a - ?' 1 * al' lr 

pr = 92 ~ hr ° 2 h,\ 2d k\ ^h r 
dt 2 r dtdr r dt r 2 ri 

ld~ , d - 2M ~ 

p = ^7^7"* + /^""r H — r^- 

/ at or r z 



The conservation (Bianchi) identity reads 

= — />' + — p + -/>' - ^P. 
The Cunningham-Pricc-Moncricf function is 



9 Dt , 9 nr , 2 „r M 



and its source term is 



2r(d~ 8~ 2- 

*odd = — ^-/it - — /i r h 

H \ or dt r 



yU \J dt Or r z 



APPENDIX D: PERTURBATION EQUATIONS IN (u,r) COORDINATES 

In this Appendix the two-dimensional Schwarzschild metric is written in the form 

ds 2 = — f du 2 — 2 dudr, 

with / = 1 — 2M/r. In these coordinates the wave operator is given by 

/ d 2 d d \ 
= -2—— + — /— *. 
\ dudr or or J 

We set A := 1(1 + 1) = /i + 2 and fj, := (I - 1)(7 + 2) = A - 2. 

The even-parity perturbations are h uu , h ur , h rr , j u , j r , K, and G. Under a gauge transformation generated by £„, 
£ r , and £ they change according to h uu — > h' uu = h uu -I- Ah uu , etc., with 



Ah uu — 2 Q £„ + £ r , 



d 2M f 2Mf ( 

' du u r 2 
_d_ _ d_ 2M_ 



The gauge-invariant fields are 
K 

The field equations are 



d 

A/i rr 2 — — £ r , 
or 

^ S Sit 7 

du 

or r 

jrt /yi ryt £i 

AG = 



„„„ _ 2 ^„ _ ™ „ + ^ + r ,|l G + „ * a _ M/ » G, 
r z r z ou 1 du Or 



Ju "r _2 ^ r ^ 



d_. _d_ 

dr ou 

Z^rjr + r 2 ^G + 2r--G 



2M 



<9 2 (9 d 

"jTjrG + r^- G-Mf G, 
auor or 



<9r J ' or 2 ' dr 
r r or ou 2 



d 2 - 23 - 19- /a- 29- \r + AM~ 
Or^ r Or r ou r Or r Or 2r 6 



Q rr 

Q u 
Q r 

Q h 



d 2 ~ 2 8 ~ r 
-K + -—K 



dudr 
du 2 



du 



<9r 



Ou 



r dr 



r > r-M d - (r-M)J g- 19- , 2/ 5_r _ / 5_r 



H — K 



2r 3 



du 



9 i 

Ou 



9 1 

Or 



9 T~, 2 l 

— K + -h ur 

or r 



dr 

2r 2 ^' 

r - M~ 



r du 



r du 



r dr 



d~ d~ d~ d~ 2{r-M)~ (r-M)f- 

h A h A K — f K - -h 4- - — h 

dr du dr r z r z 



d 2 r 



d 2 



du 2 ' rr 'dudr 
2(r -M)9- 



Q r 2 nuu 



d 2 - d 2 ~ r~M d ~ 
2^-K + f^K + ——h r 



dudr 



dr 2 



du 



roil ror 



or 



(r-M)/ fl- 2(r - M) d ~ A - 



Q» = 2h ur -fh 

r 



The conservation (Bianchi) identities read 




d O 
ou or 



9r r 2 9r 



r z r Zr 



Xr 2 - 2uMr - 4M S 
2^4 



h r 



du 



d 



dr 

ou or r 2r / 



Mf 



2M 



X 



Q uu + ^Q ur + -Q rr - —;Q r - -Q 



2r 2 



The Zerilli-Moncricf function is 
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2r 
V+2 



K + | - 2//t ur + / 2 /i rr + r^-K - rf^-K 



where A = u + 6M/r, and its source term is 

2 



24M 
A 



Q rr + 2r/Q b 
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M M 2 
- 2) + Ufa - 3)— + 84— 



fQ uu + 2Q* 



The odd-parity perturbations are h u , h r , and h^. Under a gauge transformation generated by £ they change 
according to h u — > h' u = h u + Ah u , etc., with 



Ah u 



Ah 2 



<9u ' 



9 2 
ar r 



The gauge-invariant fields are 



K 

Jl r 



h„ - 



1 d h 

2du ' 



1 d u 1 u 

K - --^-^2 + -hi- 

2 Or r 



The field equations are 



d 2 - _c^_~ _2_~ 



d 2 - 



d 



2 d - 



-P r = TT^h r - - h u + - — h u t/Iu H 5-/1 



9u 2 dudr 

d ~ d ~ J- 2M r 



<9u 



r 9u 
dr 



The conservation (Bianchi) identity reads 



= Ap« + Sipr + Ipr _ JL R 



<9u 



9r 



The Cunningham-Price-Moncrief function is 



2r ( d r 



odd 



H \dr 



fl r flu I , 



<9u 



and its source term is 



H \ ou or 



dr r 2 



APPENDIX E: PERTURBATION EQUATIONS IN (v,r) COORDINATES 

In this Appendix the two-dimensional Schwarzschild metric is written in the form 

ds 2 = - fdv 2 + 2dvdr, 
with / = 1 — 2M/r. In these coordinates the wave operator is given by 

1 / a 2 d , d \ T 

DvE,= 2 7wr + 7r/7rr- 

We set A := 1(1 + 1) = /i + 2 and M := (Z - + 2) = A - 2. 

The even- parity perturbations are h vv , h vr , h rr , j v , j r , K, and G. Under a gauge transformation generated by £„, 
£ r , and £ they change according to h vv — > = + Ah vv , etc., with 

9 3 2M 2M/ 
ov r z r z 



The gauge-invariant fields are 





— hvv 


2M 2Mf , 2 „ 
ov r z r A dv z 


-M—G-Mf—G 

av ar 


h vr 




d d 2M 2 0* 

r, Jv a Jr 2 Jr + r G + 

or ov r z ovor 


r—G + M—G, 

av ar 






a a 2 a 
or or z or 




K 


= K- 


2. 2/. a _ ^ A^ 
^ -J, - — Jr + rf—G + r—G + -G. 
r r or ov 2 





The field equations are 

Q vv = 



Q vr 
Q rr 

Q v 
Q r 



d 2 ~ 29 - 19- fd~ 2d~ \r + AM~ 
^/ r z r ar r ov r or r or 2r* 



d 2 ~ 29 - r 
X + --r-if H 



m a ~ f a ~ 



i a - 



1 r 



dvdr 

a 2 



dv 



dr 



K h rr h vv 9 h vv 

r av r ar r z 



A + 4 r 1 r 
2r^ r J 



2r 2 



. , r-M0~ (r-M)f d - 10- 2/ - /0~ 
-— # + — — # + " t-^—K + - — h vv + — — h vr - - — h vv 



dv 2 



dv 



dr 



r dv 



r dv 



r dr 



l^r + AM t _2J_~ f- ixf_~ 
2r 3 r 2 r 2 2r 2 ' 

d_~ d_~ _d_~ 2~ r- Af - 

^ tiff _ ttyf r\ f^W ~ T" o 



0W 



0r 



0r 



9- 0~ 0- 0- 2(r-M)- (r-M)f- 

/i J /i ?T _ f K A - -h A- - — h 

dv dr dv dr r A r z 



2 r 



fi 2 



2 r 



d 2 



w ^h vv + 2——K + f—K 

ar A avar 



02 

0r 2 



r-M d ~ 2d_ - _ 2 0_~ 



dv 



r dv 



r dr 



2{r-M)d~ {r-M)fd~ 2(r - M) - A- Ar 2 - 2/iMr - 4M 2 - 
- ■-■ - ti vr - - — ti rr -\- - — i\ — h, vr — — tiff) 



r 2 dr 
Q« = -2h vr -fh 

r r ' 

The conservation (Bianchi) identities read 



dr 



dr 



2r 4 



a a m 

= -^~Q VV + -rrQ vr H — tQ™ 
0t> 0r r z 



^,0^ - -Q b , 

2r 2 r 



= -^-Q vr + -^-Q rr + —J-Q vv — C —l-Q" r + — Q rr — 7r^Q r 
av ar r z r z r lr z 



- f -Q\ 



a a 
av ar 



-Q r + Q b .., 

r 2r z 



The Zcrilli-Moncricf function is 



2r 
V+2 



A 











K + -[h vv + 2fh vr + f z h rr - r—K - rf—K 



dv 



dr 
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where A = fi + 6M/r, and its source term is 
q - 4 ry 1 ni i 2 

"-"even . ^ i o\ A 

A r [jjL + 2) A 



2r 2 



dv dr 



M 



H{n - 2) + 120 - 3) h 84 



. M 2 



fQ vv - 2Q l 



fQ vv - 2Q l 



2AM 



-Q rr + 2rfQ b 



The odd-parity perturbations are h v , h r , and h 2 . Under a gauge transformation generated by £ they change 
according to h v — > h' v = h v + Ah v , etc., with 

Ah r = -— £+-£, 
or r 

Ah 2 = -2f. 



The gauge-invariant fields are 



I u 1 d h 

2 dv 

1 d 1 

K = K - -^-h 2 + -h 2 . 

2 or r 



The field equations are 



pv _ _ 



d 2 



pr 



d 2 - 29- 2 - Mr 
r aw r z r z 

2 d I , A 4 r , A*/ 1 



<9vcV <9r 2 

TT^j/ir - - - h v + --z-h v + ^h v + ^-h 

r ov 



dv 2 ' dvdr 

d~ d~ d~ 2M~ 

P = ST hr + ST h v + f^~ hr H T h r- 

ov or or r z 



The conservation (Bianchi) identity reads 



= + ^- P r + l P r _ ii p. 



dv 



dr 



The Cunningham-Price-Moncrief function is 



_2rfd~ d~ 2- 

^odd — ~o~' l v — TT-flr n v 

(j, \or ov r 



and its source term is 



/i yaw or Or r z 



r 
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